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1. Introduction 

A promising way of studying heavy quarkonium systematically in QCD is to employ an ef- 
fective field theory called potential nonrelativistic QCD (pNRQCD) ^ g]. pNRQCD is derived 
from QCD by integrating out the scale above the heavy quark mass m ^> Aqcd and the scale of 
momentum transfer mv, where v is quark velocity. 

The effective hamiltonian of pNRQCD consists of the nonrelativistic kinetic terms of a heavy 
quark and a heavy antiquark with the inter-quark potentials classified in powers of 1 /m. Up to 
0(1 /m 2 ) the effective hamiltonian has the form [^] 

H = lL + lL +V ^(r) + —V^\r) + —V^\r) 
2m\ 2ni2 «?i rri2 

+V 0) M + -V ( °' 2) to + —V^(r) + 0(l/m 3 ) , (1.1) 
m\ m z m\m2 

where m\ and ni2 denote the masses of quark and antiquark, placed at f\ and respectively. 

V^>(r=\r\), where r = r\ — T2, is the static potential. Note that ni\ and rri2 can be different. 

y(i,0)( r ) = y(o,i)( r ) (= 7(i)( r )) are the corrections at 0(l/m). W 2 '°)(r), V^- 2 \r), and V^){f) 

are the corrections at 0(1 /m ), which contain the leading order spin-dependent corrections [|], ^] 

and momentum-dependent corrections [|], Q]. 

The spin-dependent (SD) part of the 0(1/ m 2 ) correction is conventionally written as 



2m\ 2m\ J I r r J \ ^ m \ m 2 2m\m,2 

+ — f - ^> 3 (r) + f-^W , (1.2) 

m\m2 V r z 3 J 5m\m,2 



where s\ and s*2 denote the spins, and l\ = — fa = I the orbital angular momenta. Although the 
original expression of the SD correction in pNRQCD contains the matching coefficient cp it is 
assumed to be one here for simplicity. For the actual application one has to compute the matching 
coefficient in a perturbative or a nonperturbative manner depending on the matching scale between 
pNRQCD and QCD. V[(r), V^(r) are responsible for the spin-orbit corrections (fine splitting), while 
Vs(r), V$(r) are for the spin-spin corrections (hyper-fine splitting). 

The spin-independent (SI) part of the 0(\/m 2 ) correction is written as 

1 ' ■!(fiAf!P(r)}-irf ,1| W(S W6-5)+V, lUI (r)] . (1.3) 



m\m,2 \ 2 p 2r z ' 

The radial functions specified by the subscripts p 2 and I 2 , are related to the momentum-dependent 
potentials, V b (r), V c (r), V d (r) and V e (r) defined in Refs. [|, g], by 

V p(r) = V p(r) = V d (r) - \v e (r) , vf°>(r) = V^(r) = V e (r) , 

v£' 1 \r) = -V b (r) + ~V c (r), ^ U) (r) = -V c (r) . (1.4) 
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Once the functional forms of these corrections are determined in QCD, various properties of heavy 
quarkonium can be investigated systematically by solving the Schrodinger equation. 

Recently, we investigated the 0{\/m) correction [Q, and the 0{\/m 2 ) spin-dependent [[| 
[Tot l an d momentum-dependent corrections ^ in SU(3) lattice gauge theory utilizing a new method, 
and obtained remarkably clean signals up to distances of around 0.9 fm. We observed a certain 
deviation from the perturbative results at these distances. 

In this report we present our new results of the relativistic corrections at longer distances up to 



1.2 fm [11]. We then examine the scaling behavior of the corrections with respect to the change of 



lattice spacing. 

2. Formulation and Numerical Procedures 

According to pNRQCD, the 0(1 /m) and 0(1 /m 2 ) corrections can generally be expressed by 
the matrix elements and the energy gaps appeared in the spectral representation of the color-electric 
and color-magnetic field strength correlators (FSCs) on the quark-antiquark source [Jl],0|. 

Let us explain how to determine the 0(1 /m) correction as an example. 

We write the eigenstate of the pNRQCD hamiltonian at 0(m°) in the 3® 3* representation 
of SU(3) color, which corresponds to the static quark-antiquark state, as \n) = \n\f\j-ij. Then, 
the color-electric FSC, i.e. the correlator of two color-electric field strength operators E, = 
(i = 1,2, 3), attached to a quark at r\ and separated t = t\ — 1% in the time direction, takes the form 

oo 

C(r,t) = ^(OlgE^h^inlgEj^me-^^' , (2.1) 

n=\ 

where AE n o(r) = E n (r) —Eo(r) denotes the energy gap and £b(r) = (r). The 0(1 /m) correction 
is expressed with the matrix elements (0\gEi(r\)\n) and the energy gap AE n o(r) as 

vW(r) ~ - 1 & f (°\sm)\n)(n\gEj(n)\0) 

[ ) 2 d t x (AE n0 (r)f ■ iZ - Z) 

Thus, once the matrix elements and the energy gaps are known from the behavior of FSCs, one can 
compute the correction. 

We work in Euclidean space in four dimensions on a hypercubic lattice with lattice volume 
V = L 3 T and lattice spacing a, where we impose periodic boundary conditions in all directions. 
We use the Polyakov loop correlation function (PLCF, a pair of Polyakov loops P separated by a 
distance r) as the quark-antiquark source and evaluate the color-electric FSCs on the PLCF, 

C(r,t) = ((gEift^gEjft^c = ((^(n, ?1 )^(n^2))>-((^(n)))((^'(n)>) > ( 2 -3) 

using the multi-level algorithm [f7], [l(J, where the double bracket represents the ratio of expectation 
values ((•••)) = (••■)«*/ (PP*(r)}, while (••■)«» means that the color-electric field is connected to 



the Polyakov loop in a gauge invariant way. The subtracted term on the r.h.s. of Eq. (2.3) can be 



nonzero as the color-electric field is even under CP transformations. The spectral representation of 



Eq. (2.3) derived with transfer matrix theory reads [10] 



OO rji 

C(r,t) =2£(0[^(r 1 )|«)(n|g£^n)|0) e -^ r / 2 cosh((A£„ )( T -0)+O(^ (A£l0)r ) , (2.4) 
«=l l 
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Table 1: Simulation parameters used in this study. iV ts i is the number of time slices in a sublattice and A^ up( j 
the number of internal update within a sublattice, both are parameters for the multilevel algorithm [|l(J. The 
lattice spacing a is set from the Sommer scale ro = 0.5 fm. 



P=6/g 2 


a [fm] 


Ay 


spin-dependent corrections 

(L/a) 3 (7» N iupd /V con f 


spin-independent corrections 

(L/af{T/a) 7Y iupd 7V con f 


5.85 


0.123 


3 


24 4 


50000 


77 


24 4 50000 133 


6.00 


0.093 


4 


20 3 40 


7000 


33 


24 3 32 50000 100 


6.20 


0.068 


5 


24 3 30 


10000 


33 


30 3 40 50000 33 


6.30 


0.059 


6 


24 4 


6000 


39 





where the last term represents terms with exponential factors equal to or smaller than exp(— (AE i0 )T), 
which are negligible for a reasonably large T. We evaluate Eq. (2.3) via Monte Carlo simulations 
and determine the matrix element (0\gEi(r\)\n) (n\gE j(r\)\0) and the energy gap AE„o in Eq. (^J) 
by fitting, both of which are finally inserted into the definitions of the corrections, such as Eq. $L%). 
Eq. ( 2^4 ) is reduced to the form like Eq. (2T) in the infinite volume limit T — > °°. 

We define the lattice color-electric field operator, ga 2 Ej(s), from the traceless part of [t/4,-(s) — 
ulj(s)]/(2i) with two-leaf modification (an average of F^s) and F^fs — f)), where U^ v (s) is a 
plaquette variable defined on the site s. We multiply the Huntley-Michael factor Jl2| ] on the PLCF, 
Ze [p]], to the lattice color-electric field to cancel the self energies at least at 0(g 2 ). 

Other relativistic corrections can be investigated similarly. See Refs. [|| [To|] for definitions and 
further technical details. 



3. Numerical results 



We carry out simulations using the standard Wilson gauge action in SU(3) lattice gauge theory. 
Simulation parameters are summarized in Table [!]. 




> 

T3 0.5 



• P=5.85 
■ P=6.0 

▲ P=6.2 
fit curve ((3=6.0) 



r / r„ 



r/r n 



Figure 1: The 0(1 /m) correction, W(r), normalized at r = 0.8 ro (left) and its derivative with respect to 
r (right) in units of ro. The solid line is the fitting curve for the data at j3 = 6.0 with the functional form 
Vin{ r ) — —A/r 2 +Blnr + C and dV\ B (r)/dr, where the data at r/a = 2 is not taken into account in the fit. 
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Figure 2: The spin-dependent corrections V[(r), V 2 '( r )> Vj,(r) and V^(r) in units of rQ. 



3.1 The 0{\/m) correction 

In Fig. |J we show the 0(1 /m) correction, V^'(r), normalized at r = 0.8 ro together with the 
fitting curve V\ a (r) = — A/r 2 + Blnr + C. The first term is motivated by perturbation theory at 
0(a 2 ), which provides V veA (r) = — F A a ' , where Cp and Ca are the Casimir charges of the funda- 
mental and the adjoint representations, respectively, and a s = g 2 /(4n) the strong coupling [|lj]. We 
measure V ^ (r) up to r = 2.35 ro = 1 .2 fm. (r) for different j8 values show a reasonable scaling 
behavior, except for the data at r/a = 2, which are always suffered from a larger discretization error 
than the data at r/a > 3. We find that the 0(\/m) correction contains the radial behavior which 
cannot be explained by perturbation theory at 0(a 2 ). The new data set at r > 1.5 ro suggests that 



the long distance behavior is described by a logarithmic function [11, 13]. 



3.2 The 0(1 /m 2 ) spin-dependent corrections 

In Fig. ^, we present the spin-dependent corrections, V((r), V 2 '(r), Vj(r) and V^(r), where the 
data at /3 = 6.0 and 6.3 are already published in Ref. [jKj. The 0(1 /m 2 ) spin-orbit corrections, 
V[(r) and V^r), are non-vanishing up to r = 2.23 ro = 1.12 fm. The finite tail of Vjfr) is observed 
in Ref. p0|], which are again observed at further longer distances. 
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Figure 3: Test of the Gromes relation V^'(r) — V^ir) — V[(r). Comparison between the force y(°V and 
the difference of the spin-orbit corrections Vi — V[ (left), and the relative deviation from the Gromes relation 
l-(V£-V{)/VM' (right). 



The Gromes relation [ph, an important analytic relation derived from the Lorentz invariance, 

V(°)'(r)=^(r)-V/(r), (3.1) 

is approximately satisfied as shown in Fig. |[ where the comparison of the l.h.s and the r.h.s of 
Eq. (3. 1) as well as the relative deviation from the relation, 1 — (V 9 ' — V/) /V^', are plotted. If the 
Gromes relation is satisfied, this quantity should be zero at all r. We find that the deviation is 10 to 
12 % at a = 0. 123 fm, while 4 to 10 % at a = 0.093 fm, and seems to be smaller at a = 0.068 fm. 
Namely, we find a tendency such that the deviation decreases as a —> 0. 

For the spin-spin corrections, V3 (r) and V^(r), we see that they have no long range contribution. 
However a detailed analysis shows that the functional form slightly deviates from the leading-order 
perturbative expressions at intermediate distances. 

3.3 The 0(1 /m 2 ) momentum-dependent corrections 

In Fig. |J we show the momentum-dependent corrections, Vb(r), V c (r), Vd(r) and V e (r), which 
are normalized at r = 0.8 r$. Motivated by the minimal area law model [^, |J], we fit the data to the 
functional form as the static potential V(r) = —A/r+Br+C, where the data at r/a = 2 is not taken 
into account. We find that the global structure of the data are well described with this function as 
seen in Fig. |[ 

There are nonperturbative relations like the Gromes relation for the momentum-dependent 
corrections, called the BBMP relation [§, ffl], 



V b (r)+2V d (r) 



-V>(r) + - 
2 W 6 dr 



dV^{r) 



V c (r)+2V e 



rdV^jr) 
2 d~r 



(3.2) 



We examined these relations in Ref. [||] and found that they seem to be satisfied. Using the new 
data set, we are now investigating the BBMP relation carefully, which will be reported in a future 
publication. In any case, it is clear that the momentum-dependent corrections contain nonpertur- 
bative contribution as they are related to the static potential. The effect of these corrections to 
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Figure 4: The momentum-dependent corrections Vb(r), V c (r), V c /(r) and V e (r) in units of ro, normalized at 
r = 0.8r . 



the spectroscopy is not yet well examined, so that it is quite interesting to solve the Schrodinger 
equation with these corrections. 

4. Summary 

We have investigated the relativistic corrections to the static potential, the 0(1 /m) correction 
and the 0(1 /m 2 ) spin-dependent and momentum-dependent corrections in SU(3) lattice gauge the- 
ory. These corrections are important ingredients of pNRQCD for heavy quarkonium spectroscopy. 
By evaluating the color-electric and color-magnetic FSCs on the PLCF with the multilevel algo- 
rithm, and exploiting the spectral representation of FSCs, we have obtained a very clean signal for 
these corrections in the region from 0.25 fm to 1.2 fm. We have observed long-range nonperturba- 
tive contribution to these corrections, which show a reasonable scaling behavior with respect to the 
change of lattice spacing. 
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